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Aims. Focusing mirrors for X-ray telescopes in grazing incidence, introduced in the 70s, are characterized in terms of their perfor- 
mance by their imaging quality and effective area, which in turn determines their sensitivity. Even though the on-axis effective area is 
assumed in general to characterize the collecting power of an X-ray optic, the telescope capability of imaging extended X-ray sources 
is also determined by the variation in its effective area with the off-axis angle. The effective area, in general, decreases as the X-ray 
source moves off-axis, causing a loss of sensitivity in the peripheral regions of the telescope's field of view. 

Methods. The complex task of designing optics for future X-ray telescopes entails detailed computations of both imaging quality 
and effective area on- and off-axis. Because of their apparent complexity, both aspects have been, so far, treated by using ray-tracing 
routines aimed at simulating the interaction of X-ray photons with the reflecting surfaces of a given focusing system. Although this 
approach has been widely exploited and proven to be effective, it would also be attractive to regard the same problem from an ana- 
lytical viewpoint, to assess an optical design of an X-ray optical module with a simpler calculation than a ray-tracing routine. This 
would also improve the efficiency of optimization tasks when designing the X-ray optical modules. In this paper, we thereby focused 
on developing analytical solutions to compute the off-axis effective area of double-reflection X-ray mirrors. 

Results. We have developed useful analytical formulae for the off-axis effective area of a double-reflection mirror in the double cone 
approximation, requiring only an integration and the standard routines to calculate the X-ray coating reflectivity for a given incidence 
angle. The computation is easily applicable also to Wolter-I mirrors (such as those of NeXT, NuSTAR, HEXIT-SAT, IXO) and the 
approximation improves as the f-number of the mirror increases. Algebraic expressions are provided for the mirror geometric area, 
as a function of the off-axis angle. Finally, the results of the analytical computations presented here are validated by comparison with 
the corresponding predictions of a ray-tracing code. 
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1. Introduction 

X-ray telescopes have been equipped with focusing optics since 
the 70s to endow them with imaging capabilities and concentra- 
^ tion properties that enhance their sensitivity. The effective area 
k> , of the optics is one of the most important parameter determin- 
5_] ing the minimum detectable flux: to date, the X-ray telescope 
d with the largest effective area in the soft X-ray band (< 10 keV) 
is Newton-XMM with 1450 cm- per module at 1 keV, on-axis 
(Gondoin et al. |1998bT ). Beyond 10 keV, the effective area and 
the sensitivity of all focusing X-ray telescopes drop off sud- 
denly, because of the very low reflectivity of single-layer coated 
mirrors at the grazing incidence angles in use. The situation 
is expected to change with the launch of imaging hard (> 10 
keV) X-ray telescopes such as NuSTAR (Koglin et al. 120051) . 
NeXT (Takahashi et al. 2008; Ogasaka et al. l2U58] l. HEXI T-SAT 
(Pareschi i2004l l, and I XO (formerly XEUS, Parmar et al. |2004l 
Kunieda et al. I2008I I that will extend the performances of 
presently-operating X-ray telescopes beyond 10 keV, by means 
of shallow incidence angles (< 0.25 deg) and wideband multi- 
layer coatings (Joensen et al. 119951 Tawara et al. 119981 ). 

On the other side, the study of extended X-ray sources, such 
as galaxy clusters, and even cosmological surveys would benefit 
greatly from an increase in the field of view of X-ray telescopes. 
For example, eROSITA (Predehl et al. 2006 ) will provide a wide 
field of view of 61 arcmin in diameter (Friedrich et al. 120081 ) 



by using double reflection mirrors. Previous mission concepts, 
such as WF XT (C itterio et al. [T993] ) and WFI onboard EDGE 
(Piro et al. 120081 ). have been proposed to increase the imag- 
ing quality of off-axis sources, by adopting polynomial profiles 
(Burrows et al. |1992t Conconi & Campana '2 0011 ). This type 
of mirror design provides a higher resolution for off'-axis X-ray 
sources than the Wolter's, at the expense of a small degradation 
of the HEW on-axis. 

Nevertheless, not only the angular resolution, but also the 
effective area of grazing-incidence X-ray mirrors is known to 
be degraded as a source moves off-axis, due to geometrical vi- 
gnetting and the variation in the incidence angles on miiTors. 
This can have important consequences for the observation, since 
the sensitivity of the telescope over its field of view can be 
severely compromised. Moreover, for telescopes with modules 
flying in formation (as was to have been the case for SIMBOL- 
X, Pareschi et al. 120081 ). oscillations at random of the opti- 
cal module might cause an unpredictable decrease in the effec- 
tive area. Possible solutions to this problem have been studied 
(Cotroneo et al. 120091 ). although this example highlights the im- 
portance of the theoretical prediction of the effective area, on- 
and off'-axis, already at the design stage of the optical module 
development for an X-ray telescope. 

Another situation in which the theoretical computation of 
the eff'ective area is necessary occurs whenever X-ray mir- 
rors are calibrated using on-ground facilities such as PANTER 
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(Brauninger et al. 120041 Freyberg et al. I2006I I. where the X-ray 
source is located at a large, but definitely non-astronomical, dis- 
tance from the mirror. In such conditions, the measured effec- 
tive area is affected by the imperfect collimation of the incident 
beam: the correct data interpretation must account for this effect 
(see e.g., Gondoin et al. |1998b] l, on- and off-axis, to reliably re- 
construct the mirror's effective area for an astronomical X-ray 
source over all the field of view of the telescope. 

While the prediction of the on-axis effective area for a 
double-reflection, grazing-incidence mirror is rather simple, it 
becomes difficult in general for a source off-axis. Such evalu- 
ations have been completed so far using ray-tracing codes (see 
e.g., Mangus & Underwood 1969; Zhao et al. ,2004) , starting 
from the optical design of mirrors and the structure of the re- 
flective coating. These codes are in general very efficient, but 
time-consuming, computationally intensive, and affected by sta- 
tistical errors related to the number of rays that can be traced. It 
would therefore be beneficial to develop a method for computing 
the efifective area analytically, which would ease the assessment 
of the effective area and, among other things, the optical design 
of a wide field X-ray telescope. 

Few suitable analytical tools have been developed to date. 
Van Speybroeck and Chase (|1972| l found - by means of ray- 
tracing - that the geometric, collecting area of a Wolter-I mirror 
decreases with the off-axis angle 6 of the source, with respect to 
the on-axis geometric area Aoo(O), according to 

Aeo(0)=Aoo(O)(l-^), (1) 

where oq is the incidence angle for a source on-axis. However, 
we lack a general, analytical method to derive the effective area 
of an X-ray mirror with a given reflective coating, as a function 
of the off-axis angle of the X-ray source. 

In this paper, we present a solution to that problem. We de- 
velop an analytical approach that can be applied to double cone 
grazing-incidence X-ray mirrors and, with reasonable accuracy, 
to Wolter-I mirrors (unless the f-number is small). The limits of 
this approximation are discussed in Sect. |2l In Sect. [3l we de- 
rive general integral formulae (such as Eq. (|25] |) to compute the 
off-axis effective area for a double-reflection X-ray mirror with 
shallow incident angles, for any reflective coating. As a partic- 
ular case, in Sect. |4] we obtain some algebraic expressions for 
the geometric area and verify that the well-known Eq. ([TJ can 
be derived as a particular case. In Sect. |5] the predictions of the 
analytical approach are validated, for some particular cases, by 
means of a comparison with the outputs of a ray-tracing routine. 
The results are briefly discussed in Sect. |6l 

We note that we assume that the off-axis mutual obstruction 
of mirrors in densely nested mirror assemblies has a negligible 
effect. Therefore, the results are valid for either isolated dou- 
ble cone or Wolter-I mirrors, or for mirror modules known to be 
negligibly obstructed, such that their effective area simply equals 
the sum of the contributions of the individual mirrors. The quan- 
tification of the off-axis obstruction in mirror assemblies will be 
considered in future. 

2. The double cone approximation in the 
computation of the effective area of a Wolter-I 
mirror 

We consider, in a preliminary way, a grazing-incidence Wolter-I 
mirror, and an on-axis photon source (Fig. (TJ. The optical axis 
is aligned with the z axis. We define Rm to be the radius at the 




Fig. 1. Meridional section of a Wolter-I mirror, here in the case 
L2 > Li. a) source on-axis, at astronomical distance; b) source 
still on-axis, but located at finite distance (all angles are greatly 
exaggerated). 

parabolic end (i.e., the maximum radius), Rq the radius at z = 
(the intersection plane), R„, the radius at the hyperbolic end (i.e., 
the minimum radius), F the focal point, and / the distance of 
F from the intersection plane (i.e., the focal distance). In gen- 
eral, we refer to "primary" and "secondary" segments, instead 
of "parabola" and "hyperbola". We denote with Li the primary 
segment length along the z axis, and L2 that of the secondary. 
The polar coordinate is (f. 

Because of the surface curvature, the incidence angles on 
the two surfaces vary in general with the z coordinate. We de- 
fine a{z) to be this angle on the primary (0 < z < Li) mir- 
ror segment, and on the secondary (-L2 < z < 0) mirror 
for a source on the optical axis, at infinite distance. In general 
(Van Speybroeck and Chase ,1972) , the incidence angle of rays 
close to z = 0, ao, is the same for both surfaces. Therefore, we 
have the well-known relation 

/Jo =/tan(4ffo). (2) 

Since ao is shallow, we hereafter assume that the tan(4Q'o) func- 
tion can be approximated by Aao itself. 

To obtain an analytical expression for the effective area of a 
Wolter-I mirror when a source is off-axis by an angle 6, we have 
to find algebraic expressions for: 

1. the variation in the incidence angles, on which the mirror 
reflectivities depend, over the primary and secondary mirror 
segment surfaces; 

2. the collecting area of the primary segment; 

3. the fraction of rays reflected by the primary segment that are 
incident on the secondary. 

To compute the effective area for a source off-axis, we begin with 
the prototypical case of a mirror with a source on-axis, but at a 
finite, although very large, distance D, which is the usual con- 
figuration for on-ground calibration facilities. All rays impinge 
the primary segment within a meridional plane, and in the dou- 
ble cone approximation they are all incident at the same grazing 



D. Spiga et al.: Analytical computation of the ofF-axis Effective Area 



3 



angle, independently of z and ip. The finite distance of the source 
causes the beam to have a divergence at the primary segment, 
which can be assumed to be constant as long as D » Li, with 
an half-aperture angle of 6 - Rq/D. In this simplified case, ev- 
ery point of the mirror sees the source off-axis by an angle 6, 
regardless of cp. We demonstrate hereafter that the solution to 
these problems is simpler to express analytically if the profile of 
the mirror can be approximated with a double cone, if we are 
interested only in the effective area. In contrast, the curvature of 
mirrors along the axis is essential for their angular resolution, 
but we do not consider this aspect here. 

In the remainder of this section, we quantify the errors 
caused by the substitution of a Wolter-I profile with a double 
cone, by keeping Rq, ao, Li, and L2 constant. The problems to 
be faced are related to the system geometry, rather than to the ab- 
solute size of the mirrors. Hence, it is convenient to report the re- 
sults of this section in terms of both the f-number, /# = f/{2Ro), 
instead of / itself, and L'j = Li /(27?o), LL = L2/{2Ro), which are 
normalized mirror lengths. Using Eq. /# can be written as 



8ao 



(3) 



We mention that the true focal length of a double cone is slightly 
longer than that of a Wolter profile with the same Rq and ao, 
because of the different focusing properties. In all cases, we al- 
ways refer to / as the focal length of the corresponding Wolter-I 
mirror, to ensure that Eq. (|2| retains its validity. 

2.1. The incidence angles 

For a double cone profile, the incidence angles for an on-axis 
source at infinity are constant along a meridional plane, and 
equal to ao. For a Wolter-I profile, the slope of the mirror pro- 
file changes slowly with the z coordinate. Because of the con- 
cavity of profiles, a(z) < ao on the primary segment, and 
a(z) > ao on the secondary segment. We assume for simplic- 
ity that Li - L2 - L, and define the (positive) angle variations 
Aai ^ ao - a(+L) and Aa2 - a{-L) - ao. For a Wolter-I profile, 
the ratio AQ'2/Aci'i varies between 1 and 2, going from large to 
short /#. For this estimation, we assume that Aa2 ~ Aai and we 
simply denote their value with Aa. 

Now consider a ray from an on-axis source at infinity 
(Fig. [T^). After being reflected firstly at z = +L, it undergoes 
the second reflection at z ^ -L, at a radial coordinate 7?,„, and 
is focused at z = -/. The total photon deflection at z = -L is 
2(0-0 + Ao-i) + 2(ao + Aq;2); therefore. 



R,„ = 4(/ - L)(ao + Aa). 
For small ao we also have that 
R„, - Ro - 3Lffo- 



(4) 



(5) 



By comparing the last two equations, using Eqs. (|2|i and (O, with 



the definition of L' 



L2 



L' , we derive the maximum variation 



in the incidence angles over the mirror profiles, 
Aa ^ U 

~^ ^ Hf#-L'y 



(6) 



This equation provides the error in the incidence angles intro- 
duced by the double cone approximation. Such an error is com- 
puted for some practical cases and reported in Table \T\ For ex- 
ample, we note that for NeXT-HXT and SIMBOL-X the error is 
only 0.4%. This means that, assuming ffo - 0.27 deg at most. 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

2 3 4 5 6 7 8 9 10 11 12 13 
f-number 

Fig. 2. Percent accuracy of the primary segment cross-section 
area, as seen by a source on-axis from infinity, in approximating 
a Wolter-I profile with a double cone. 

Table 1. Some examples of approximation introduced 
by the double cone geometry (optical parameters after 
Pareschi et al. 20()8] Ogasaka et al. I2b08l Koglin et al. 120051 
Gondoin et al. 1998al l. 





SIMBOL-X 


NeXT-HXT 


NuSTAR 


XMM 


/ 

L 


20 m 
324 mm 
300 mm 


12 m 
225 mm 
200 mm 


~10m 
~169 mm 
200 mm 


7.5 m 
346 mm 
300 mm 


Aa 

ao 


0.4% 


04 % 


0.5% 


1% 


AAi 
Al 


-0.18% 


-0.20% 


-0.25% 


-0.5% 


^(6 = 0) 


1.5% 


1.7% 


2% 


4% 



the angular variation along the profile is Aa ^ 4 arcsec, which 
has almost no practical influence on the reflectivity of wideband 
multilayers. 

We now consider the source on-axis at finite distance, with 
6 < ao- For a double cone (Fig. [TJ)), the incidence angle on 
the primary segment becomes ai = ao + for all if. Similarly, 
the incidence angle on the secondary segment is also constant 
and equal to a2 - ao - 6, provided that ao > S. Within the 
approximation of Eq. (|6|l, this is also true for a Wolter-I mirror. 



2.2. The maximum diameter 

In terms of effective area, another concern of approximating a 
Wolter-I profile with a double cone of the same Ro, ao, Li, and 
L2, is the different cross-section of the primary segment. For a 
double cone profile, the maximum diameter is Rm = Rq + aoL\ . 
For a Wolter-I mirror, because of the profile concavity, the max- 
imum diameter is smaller by an amount that depends on L' and 
/#. Consequently, the geometric area of the primary segment is 
smaller for a Wolter profile than for a double cone with the same 
size and ao. However, the (negative) error caused by the substi- 
tution of a Wolter-I profile with a double cone, defined as 



AA I _ A I (Wolter) - A 1 (double cone) 
"a^ ~' A 1 (Wolter) ' 



(7) 
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8 = 1 arcmin 




10 -1 



L' = 0.5 



Fig. 3. Vignetting due to the finite distance of the source on-axis. 



as a function of the f-number, for L'. 



L'^ and 4 diff'erent val- 



ues of 6. Solid lines are computed assuming a Wolter-I profile, 
dashed lines a double cone (Eq. (fTTTi). The results for the double 
cone are completely independent of the choice of L': those for 
the Wolter-I are computed assuming L' - 0.5. 

is, in general, of a few percent as far as L' < 1 (see Fig. |2]i. 
In addition, its magnitude decreases rapidly for increasing /#, 
since it is expressed well by the empirical formula 



(8) 



AAi L' 

of the same order of magnitude as the expression in Eq. (|6]). 
This approximation is also computed for some telescopes and 
reported in Table [T] 

For an on-axis source at a finite distance, the area of the cir- 
cular corona between 7?^ and R()-Li6 should be considered: the 
area is thereby increased by the same term ~ 2nRf,L\6, for both 
Wolter and double cone. Hence, the approximation of the area is 
even better than the expression in Eq. ([8]). 

2.3. Geometric vignetting for double reflection 

We now consider the third problem. The case of a source on-axis, 
at a finite distance, represents the simplest case in which there is 
geometrical vignetting for double reflection. If the source is on- 
axis, at infinite distance (with L2 > Li), all rays undergo two 
reflections and reach the focus of the mirror. If the distance is fi- 
nite, not all photons reflected by the primary segment also hit the 
secondary: for instance, photons reflected at z ~ +Li can miss 
being reflected by the secondary, so they are prevented from be- 
ing focused, as depicted in Fig.[T] We can quantify this geometric 
effect by means of the vignetting factor V, which is defined to be 
the ratio of the number of photons reflected twice to the number 
of photons reflected by the primary segment (assuming, ideally, 
that the mirror reflectivity is 1). In this configuration, if the mir- 
ror profile is a double cone, we can compute V using the simple 
formula (see AppendixlAlfor a derivation) 



V ^ 



L] at 



(9) 



valid if L2a2 < iici, otherwise V = 1. If either cci < or 
a2 < 0, then V = 0. In Sect. |3] we see that this result can even 
be generalized to the case of a source off-axis. 



V 6 = 1 arcmin 

O 8 = 3 arcmin 

□ 8 = 5 arcmin 

A 8 = 10 arcmin 




' I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

5 10 15 20 25 30 35 40 45 50 
f-number 

Fig. 4. Percent accuracy of the double cone approximation in the 
estimation of the vignetting factor for a source on-axis at finite 
distance, as derived from the curves of Fig. [3] (symbols). Lines 
are traced using the empirical Eq. (fTSli with L' = 0.5. 



If we consider that in this case ai - ao + 6 and a2 - ao- 6, 
an alternative expression for Eq. (|9]l is 



V 



L2(D-4f) 
Li {D + AfY 



(10) 



where we have used the definitions of oq and S. Using Eq. ([3]l 
and normalizing the L's to the mirror diameter, we can also write 
Eq. © as 



V 



L'^(l-Uf#) 
L\{l+%5f#y 



(11) 



It is interesting to note that a formula similar to Eq. ( fTOb for 
the effective area was empirically found by Van Speybroeck and 
Chase ( |1972| l. but with a 2.5/ term instead of 4/. 

Rigorously, Eqs. (|9|l to ( fTTT i hold only for a double cone mir- 
ror, but they are expected to be applicable also to Wolter-I mir- 
rors of a sufficiently large /#. To determine the approximation 
that can be achieved, we computed by ray-tracing the exact V 
factor of a Wolter-I mirror as a function of /#, in the case of 
L'j - L'^- 0.5, for fixed 6 values. In Fig.[3l the results are com- 
pared with the findings of Eq. ( fTTT i. The vignetting for the double 
cone is independent of L' , whereas for the Wolter's it depends 
markedly on L' only for small /#. 

Inspection of Fig. [3]shows that 

- for a fixed 6 and not too large /#, V is larger for a Wolter-I 
mirror than for the corresponding double cone; therefore, the 
double cone approximation returns, in general, pessimistic 
predictions of the eff'ective area. 

- In double cone approximation, V < 1 for all /#, i.e., there 
are always lost photons for double reflection. With a Wolter- 
I profile, y is 1 when /# becomes sufficiently small. 

- For a fixed 6, the predictions of Eq. (|9]l approach the ex- 
act calculations as /# is increased. For large /# and large S, 
where V ^ 0, the V factor of Wolter becomes slightly smaller 
than that of the double cone. 

From the curves of Fig. |3] we computed the error in the 
vignetting factor when we apply Eq. (|9l) to a Wolter-I profile, as 
in the definition of Eq. ([8]l. This quantity is plotted in Fig.|4]for 
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the case L' - 0.5. It can be seen that, in the cases interesting for 
us, i.e., for sufficiently large /# and not too large values of 6, ^ 
is positive and of the order of a few percent. We repeated this 
exercise for several values of L' in the interval 0.25 - 1 to in- 
vestigate the dependence of ^ on this parameter, and it turned 
out that this ratio, if positive and if y(Wolter) < 1, can be ap- 
proximated very well by the empirical formula (see Fig. H] and 
Table [B 

AV , / 1 \ 

-y-L'[j-,-y^y (12) 

where y ^ 14.3 is with very good approximation a constant in 
the explored range of L' values. This equation has the same kind 
of dependence as Eq. (|6]l for the slope variation along the profile 
of the mirror and Eq. (O for the area of the primary segment. 

For small /# values, Eq. ( fT2] ) is not obeyed because of the 
saturation of y(Wolter) to 1, so the expression of ^ should be 
interpreted as an upper limit. Finally, for very large 6 the error 
deviates from Eq. (fTST i as it becomes negative, but this occurs 
only when V — > 0, so its weight in determining the effective area 
is expected to be negligible. 

We can therefore conclude that the error introduced by 
the double cone approximation, regarding the effective area, is 
definitively smaller than L' /f#. In other words, the double cone 
approximation is valid when L «: /, a condition fulfilled in al- 
most all practical cases. Within the limits of this approximation, 
we derive in the next section the effective area for a source on- 
and off-axis. 



3. The on- and off-axis effective area of a Wolter-I 
mirror 

3.1. On-axis source 

We assume that we can approximate the Wolter-I profile with 
a double cone profile, by adopting the tolerances estimated in 
Sect.|2] In the following, we adopt the convention to denote with 
AoiA, 0) the effective area of the mirror at the photon wavelength 
A, for a source off-axis by 0, at a distance D (finite or infinite). 
When referring to the geometric area, we omit A and use the 
notation A/)(0). Firstly, we assume L\ - L2 - L. The geometri- 
cal, collecting area of the primary segment, as seen by a source 
on-axis at "infinite" (i.e., astronomical) distance is 

Ai,<^(0) = n{Rl^-Rl) - 27tRo(Rm-Ro) - 2nRoLao, (13) 

where we use the approximation Rm - Ro - La^. In this case, 
all reflected rays also undergo an identical reflection on the sec- 
ondary segment (Fig. [1), therefore Aoo(O) = Ai k,(0). We now 
denote by r^{a) the reflectivity of the mirror at the photon wave- 
length A, for a generic incidence angle a. The form of this func- 
tion depends on the coating structure: for a single layer coating, 
which operates in total external reflection, it slowly decreases 
up to the critical angle for A, followed by a sudden cutoff. If a 
multilayer coating is used, r^{a) is a more complicated function 
and can be computed using one of the standard methods (e.g., 
Parrat |T954l Abeles[T950jl, by also including the effect of rough- 
ness using, e.g., the Nevot-Croce (1 19801 1 approach. 

Since photons are reflected twice at the same angle ao, we 
multiply the geometrical area by the squared reflectivity to ob- 
tain the effective area at the photon wavelength A 

A^{A,0) = Ae<,(0) • rjiao) = InRoLao ■ r^iao). (14) 
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This is a well known result. We now keep the source on-axis, but 
at a finite distance D and assume more generally that Li L2 
(with Li = L2 as a particular case). As already discussed in 
Sect.|2] all mirror sectors see the source off-axis by the same an- 
gle 6 - Ro/D. The effective area of the primary segment thereby 
becomes 

Ai,o(i, 0) = InRoUai ■ r^(Qfi), (15) 

where a\ =0-0 + <5- The effective area for a source at finite dis- 
tance is then obtained from Eq. (fTSl l. times the vignetting factor 
of Eq. (|9]l, times r^{a2), the reflectivity of the secondary mirror 
segment, 

Ad{A,Q) = InRoLiaiV- r,i(a,)rAa2), (16) 

where a2-ao-6.IfV<\, substitution of Eq. (|9]l yields 

Ao(/l,0) = ItiRo ■ L2a2 ■ r^(ai) r,i(Q'2), if L2a2 < Li ai. (17) 

As one might expect, the geometric area is that of the secondary 
segment, projected onto the wavefront after the primary reflec- 
tion. We note that by setting <5 = we retrieve the on-axis result, 
in Eq. ( fT4l i. On the other hand, if L20'2 > L\a\, all the primary 
segment is effective in the double reflection, so V = 1 . This may 
occur with a divergent source on-axis, if L2 » Li . In the absence 
of a geometrical vignetting, the effective area becomes 

Ad{A,Q) = InRo ■ Liai ■ r^iai) r,^{a2) if L2a2> Uai. (18) 

Comparison of Eqs. (fTTj l and ( fTSl l indicates that the effective area 
can be written as 

Ad{A,0) = InRo mm{Liax, 120:2) r;^{a{) rA{a2), (19) 

which represents the general expression for the effective area 
seen by a source on-axis. 

3.2. Off-axis source: integral formuia 

We can now compute the effective area for a source off-axis. 
Because of the axial symmetry of the mirror, we choose the x 
axis so that the source lies in the xz plane, on the side of the 
positive X axis (refer to Fig.|5j. We define > to be the angle 
between z and the source direction. 

With respect to the on-axis case, there are some additional 
difficulties. The ray no longer lies within a meridional plane of 
the mirror, so the polar angles of the impact positions on the pri- 
mary and secondary mirror segment, tpi and (p2, differ in general, 
and the incidence angles ai and a2 also vary with them. 

Nevertheless, since the maximum distance of the two impact 
points is ~ Li -H L2 ~ 2L, the off-plane linear displacement is 
26L at most: therefore, \ifi2 -ifi\< AOL', which is in general 
negligible with respect to ip\ and 1^2 themselves. This result can 
also be derived more rigorously (see Appendix iBli. 

In the limit of small 9, we can henceforth assume with good 
approximation (f\ x (f^: for this reason, in the remainder of this 
section we denote with ip the nearly-common value of the two 
angles. In Appendix iBl we also derived other important results: 

1 . analytical expressions for ai (ip) and a2{'f) can be found with 
some algebra. For a small angle 6, these functions reduce to 
the simple expressions 

ai{if) ^ ao + S-6cos(p (20) 
0-2(1^) - ao - S + Ocostp, (21) 

which generalize the expressions for the incidence angles of 
Sect. 13. II to the case of an off-axis source; 
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to the focal plane 



Fig. 5. A grazing-incidence double cone mirror, illuminated by 
an ofF-axis source: the dashed line is a ray path. The optical axis 
lies along the z direction and the x axis is chosen for the source 
5, at a distance D - Ro/6, to lie in the xz plane. The source di- 
rection forms an angle with the z axis. The azimuthal positions 
of the reflection points are located along with the tpi and (p2 polar 
angles. 



2. always in the limit of small ao, 6, 9, the vignetting for double 
reflection (Sect. l23b can be calculated using a generalization 
ofEq. I®, 



(22) 



Here the V factor varies with ip because it depends on a\ 
(Eq. (Unil) and 02 (Eq. dlB)- As in Eq. Eq. ^ is valid 
only if < V{ip) < 1. If V exceeds 1 at some ip*, then 
all the primary segment sector is effective and Visp*) = 1 . If 
either ai or ao is negative at some (f*, no reflection occurs 
on the optical side of the mirror and Visp*) - 0. 

We now draw our attention to a small mirror sector between 
<p and ip+A(p, when A(p — > 0. Equations ( l20t and (ISTT i can then be 
used to compute the local incidence angles, and Eq. ( |22] | to com- 
pute the local vignetting factor. We assume, initially, that ai > 
and a2 > 0; therefore V > 0. If V((p) < 1, we can repeat the pas- 
sages in Sect. 13. II for a mirror sector at (p according to Eq. ( fTTb . 
Otherwise V((p) - 1, all the primary segment sector is effective, 
so we can write the sector area according to Eq. (fTSl l. Finally, we 
obtain for the effective area of the infinitesimal sector. 



AAd(A, 0, if) = Ro (La)^in r^(ai) ri(ff2) A<p, 



(23) 



where (La)ini„ - min[Liai((p), L2Q'2(v')]- InEq. ( |23] ), we omitted 
the explicit dependence of a[ and 02 on (p: in the following, we 
adopt the same convention to simplify the notation. 

The expression for (La)^i„ is correct as long as the incidence 
angles are non-negative. If either ai or a2 become negative at 
some (p, no contribution to the effective area can be given by that 
sector, so (LQ')min is zero at that ip. We can then write without 
restrictions that 



(La), 



max (0, mm[Liai((p), L2a2((p)]} ■ 



(24) 



The effective area of the entire mirror is then given by integration 
over (p 



h r (La) 
Jo 



Ad(A,0) = 2Ro (La)mmrA(ai)rA(a2)d(p, 
Jo 



(25) 



where the factor of 2 comes from the symmetry of the mirror 
with respect to the off-axis plane. In the particular case Li = L2, 
Eq. (|25] | becomes 



Ad(A. 



,0)^2RoL [ a„ 
Jo 



rA(ai)rAa2)dip, 



(26) 



where a^ain - ^^n[ai((p),a2((p)] if positive, and zero otherwise. 
Using Eqs. (|20] | and i2T[ . we can conveniently write o-min in the 
compact form 



max(0, Q-Q - |i5 - 6 cos ip\). 



(27) 



which is valid also for 6 > ao. Finally, for a source at infinity 
there is additional symmetry with respect to the y axis, and, if 
< ao, we can write Eq. ( |26] l in an even simpler form. 



•I 



■71/2 

Aoo(A, 0) = 4-RoL I (ao - 0cos ip) o(ai) o(q'2) d^- 



(28) 



4. Applications to the geometric area 



An upper limit to the effective area Ar,(A,0) is represented by 
the geometric area, Ao(0), obtained by simply setting the mir- 
ror reflectivity to 1, for all A and a. In the most common case, 
L\ - L2, the expression of the geometric area is 



Ad(0) 



- 2RoL I a„ 
Jo 



(29) 



where amin is given by Eq. (l27l l. In this section we solve this 
integral explicitly and provide analytical expressions for Ab(6'). 

4.1. Source at infinity 

For a source at infinity with < ao, it is more convenient to use 
Eq. (ESl) with rA(a) = 1 



•I 



■k/2 

Aoo(0) - 4-RoL I (ao - 0cos <p) dip. 



(30) 



which can be immediately solved 

AU0) = 2nRoLao - 4-RoL0. (31) 
Recalling Eq. ST4i for the on-axis area Aco(O), we obtain 



0<an i 20 

A^{0) =°A,o(0) 1 

nao 



(32) 



which is exactly Eq. ([T]i found by Van Speybroeck and 
Chase (I1972I I. after approximating tt ^ 3. 
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For off-axis angles larger than oq, Eq. ([32]) is no longer valid. 
An extension of the curve A^oiff) for > uq can be obtained from 
Eq. ( l30t , after setting the integrand to when cosi^ > q-q/^, 
according to Eq. ( |27] l. The result is a non-linear function of the 
6/ao ratio. 



(0) 







n 





02 ffO 

- 1 + arccos — 

^2 o 







(33) 



which is identical to Eq. (DZb for 9 - oro; this is correct because 
the geometrical vignetting must be a continuous function of 6. 
We note that, for sufficiently large 6, Eq. (l33T l can be approxi- 
mated well by Aco(O) ao/ (nO). 

A complete vignetting curve for 5 = is shown in Fig. |6] 
(solid line). The predicted deviation from linearity is verified in 
Sect.|5]by means of an accurate ray-tracing routine. 

4.2. Source at finite distance 

We now consider the variation in the geometric area, for a source 
at a finite distance (6 > 0). In this case, the integration depends 
on whether 6 < ffo/2 or not. 

We consider firstly the case 6 < ao/2. We then assume ini- 
tially 6 < S. With these conditions, 6 cos tp < 6 < ao - 6 for all ip, 
therefore Eq. dZTl i becomes 



a'mm(v) ^ ao - 6 + 0COS(p. 



(34) 



Substituting this expression into Eq. ( l29b and solving, we derive 
the area 



AD(0) = 2nRoL(ao-d), 



(35) 



that is, Ao{0) - Ad{Q) (see Eq. (fTTb). In other words, the mirror 
geometric area is constant as far as 6* < 5 < ao/2. This is often 
observed in optics calibrations at on-ground facilities (see, e.g., 
Gondoinet al. 1998b). 

We now increase beyond 6. Since 6 < ao-6hy hypothesis, 
we can consider the case 6 < < ao - 6 < ao + 6. We are 
therefore allowed to write Eq. ( |29] ) as 



y^arccos | 

ioiO) = 2RoL I (ao + 6 - cos (f) dip ■ 
Jo 

+ 2RoL I (ao-S + cos ip) dip, 

«7 arccos ^ 



where the integrands are always positive. This yields 
Ad{0) = A„(0) 



1 



25 
nao 



arcsin - + 



02 
52 



1 



(36) 

(37) 
and for 



We note that for 5^0, Eq. (|37] ) reduces to Eq 
0^510 Eq. ( I35]) . as expected. 

We now suppose ao - 6 < < ao + 6 . In this case, the 
integration returns some more terms 



Ad{0) = AU0){1- — 
nao 



arcsin - + 



52 



1 

nao 



V6'2 - (ao - 5)2 - (ao — S) arccos — 



(38) 



This equation, as expected, returns the same result as Eq. (|37T ) 
at = ao - 5, and is valid for ofF-axis angles in the interval 



1.0 -1 
0.9 - 
0.8 
0.7 

g 

o 

< 0.4 
0.3 - 



8 = deg 
8 = 0.04 deg 
8 = 0.08 deg 
8 = 0.12 deg 
8 = 0.1 6 deg 
8 = 0.22 c 



0.0 



0.1 



f""" I i i ' I ' I ' I 

0.00 0.04 0.08 0.12 0.16 0.20 
Off-axis angle (deg) 
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0.28 



Fig. 6. Normalized geometric area, Ad(0)/Aoo(Q), for a double 
cone mirror with a q-q = 0.2 deg, as a function of the ofF-axis an- 
gle, for different S angles due to the finite distance of the source. 
The curves are traced using the analytical formulae reported in 
Sect.m 

[ao - 6,ao + S\. Therefore, for (5 — > it reduces to Eq. ( [32] ) 
and ( [331 for the single point - ao. 

Now, consider the case 6 > ao/2, unlike we have hitherto 
assumed, but still that 6 < ao. This time ao - 5 < 6: the condi- 
tion < 6 is insufficient for avoiding negative values of o-min in 
Eq. ( [34] ), therefore Eq. ( [35] ) is valid only if < ao - 6. Beyond 
this limit and up to = 6, Eq. ( [37] ) should be replaced by 



Ad(0) = 



Ad(0) 



6-ao 



02 



^ \ (ao - S)^ 



- 1 



(39) 



which, unexpectedly, is an increasing function of in the interval 
ao - 5 < < 6. When 6 < < 6 + ao, the geometric area again 
follows Eq. ([38] ). which for 6 > ao/2 exhibits a maximum at 



^2 + "^aod - ^al : 



(40) 



such a maximum is not present if 6 < ao/2. 

In a similar fashion, one can easily compute the total geo- 
metric vignetting for > ao + 6, even if the resulting expression 
would be too long to report here, and it is also possible to de- 
rive similar expressions for the unusual case (5 > ao. To provide 
the reader with a qualitative visualization of the overall trend of 
the expressions we just derived, we traced in Fig.|6]some curves 
of total vignetting as a function of 0, for ao = 0.2 deg and dif- 
ferent values of 6. To do this, we used Eqs. ( [32] ) and ([33] ) and 
Eqs. ( [35] ) to ( [39] ) in the respective intervals of validity, and also 
the expression for > ao + S, which is not reported here. 

To summarize the results presented in this section, from 
Fig.|6]we can draw some conclusions: 

- if 5 = 0, the geometric area decreases linearly up to = ao, 
then approaches zero as is increased beyond ao; 

- if < 5 < ao/2, the geometric area remains constant for 
< 6, then decreases monotonically to zero; 

- if ao/2 < 6 < ao, the geometric area remains constant for 
< ao- 6, then increases, reaches a maximum and eventu- 
ally decreases to zero; 
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symbols: ray-tracing 



arcmin 



Ray-tracing: 



— I — ' — I — ' — I — ' — I — ' — i— 

0.00 0.04 0.08 0.12 0.16 0.20 
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Fig. 7. Comparison between some analytical curves of Fig. |6] 
(lines), and the results of an accurate ray-tracing (symbols) for a 
Wolter-I mirror with the same ao (0.2 deg) and L' - 0.5. To avoid 
confusion in the figure, not all curves were represented. Note the 
very good matching of the two methods, within the error bars of 
the ray-tracing. 

- if 6 > ao, the geometric area is zero for < 6 - oq, then 
increases, reaches a maximum and eventually decreases to 
zero. 



5. Validation with ray-tracing results 

We verify the formulae derived in previous sections by means 
of a comparison with the results of a ray-tracing routine. To 
make the comparison easier, we consider the case of a single 
Wolter-I mirror. A first point to be checked is the geometric area, 
which was extensively analyzed from the analytical viewpoint 
in Sect. |4] In particular, the curves of Fig. |6]are easily verifiable 
along with a ray-tracing routine, without any assumptions about 
the reflective coating. The comparison is interesting especially 
when testing the analytical formulae in the non-linear regions, 
which are usually not exploited in X-ray optics. 

We displayed in Fig.|2]the comparison between 3 representa- 
tive analytical curves of normalized geometric area (from Fig.|6]l 
and the findings of a ray-tracing run on a Wolter-I mirror with the 
same value of oq (0.2 deg), corresponding to /# « 36. For the 
ray-tracing simulation, we adopted the reasonable value L' - 0.5, 
whereas the analytical curves are largely independent of this pa- 
rameter, as long as the double cone approximation is applicable. 
More exactly, from the discussion of Sect. |2] the approximation 
that we introduce is not larger than L' /f# = 1.5%. As can be 
noted, the agreement between the two methods is complete, also 
for large off-axis angles where the vignetting ceases to be linear, 
and even for 6 > a^. On the other hand, the analytical computa- 
tion, with a nearly-continuous 6 sampling, corresponded to a few 
lines of IDL code, whilst the ray-tracing routine is a complex 
program, which requires the acquisition of more than lO'* pho- 
tons to reach a statistical error of a few percent for each point. 

We now wish to compare the predictions of Eq. ( |26] l for the 
effective area, as a function of the photon energy, with the re- 
sults of the ray-tracing. As a first case, we check the results 
for the largest mirror shell of a Newton-XMM optical mod- 
ule (Gondoin et al. [T998b), with / = 7.5 m, Rq = 346.2 mm 
(ao - 0.66 deg), and L - 300 mm. The errors in the primary mir- 
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Fig. 8. Comparison of effective areas of the largest Wolter-I mir- 
ror of Newton-XMM, as computed from ray-tracing (symbols) 
and from Eq. ( l26b (lines). The source is supposed to be at infinity 
(6 = 0). The agreement is very good to within a few percent. 



ror segment area and the incidence angle variation along the pro- 
file introduced by the double cone approximation are less than 
1%, and the error in the vignetting for double reflection is 4% at 
most (see Table [T]), so Eq. (|26] | can be reliably applied to com- 
pute the effective area of the mirror, both on- and off-axis. 

The comparison between the results obtained from a ray- 
fi-acing (symbols) and from the application of Eq. ( |26] l (lines) 
is shown in Fig. [8] for a source at infinity and 4 different val- 
ues of the off-axis angle. For both computational methods, the 
reflectivity of the Gold coating was computed by assuming the 
same surface roughness rms of cr = 4 A, regardless of the photon 
energy and the incidence angle. This is not completely correct, 
because the spectral window of the roughness power spectrum 
effective for X-ray scattering changes with A, a^, 02, and the 
size of the region over which the image is integrated. A vari- 
able cr, computed from the power spectrum of roughness in a 
variable frequency range, should be adopted (Spiga et al. |2009l l. 
Nevertheless, this has no relevance to the present comparison 
and we retain cr as a constant. 

As can be noted from Fig. [8] the comparison provides a 
good agreement between the two methods, at all considered 
energies and off-axis angles. The analytical method underesti- 
mates the ray-tracing findings by only 3.7% at most, as fore- 
seen, an amount close to the statistical error for the ray-tracing. 
This confirms the correctness of the analytical formula (Eq.(|26ll) 
for Wolter-I mirrors within the approximation limits stated in 
Sect.m 

Finally, as an application to the hard X-ray band (> 10 keV), 
we consider the case of a long focal length (/ - 20 m) X-ray mir- 
ror with a wideband multilayer coating. Because of the complex 
dependence of multilayer reflectivity on the photon energy and 
the incidence angles, this example places the analytical method 
to the test, because any departure of the incidence angles from 
the true values would result in a displacement of the reflectance 
peaks. As simulation parameters, we assumed - 296.2 mm, 
and L\ - L2 - 300 mm, ckq = 0.106 deg. The geometrical area 
for a source on-axis, at infinity, would be 5.17 cm^. The ratio 
expressing the departure of the double cone from the Wolter's, 
L' I f#, is only 1.5%. The multilayer coating is supposed to con- 
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Fig. 9. Comparison of effective areas of a multilayer-coated 
Wolter-I mirror, as computed from a ray-tracing (symbols) and 
from Eq. (|26] | (lines). The source is supposed to be at infinity 
(6 = 0). The error bars of the ray-tracing outputs are not shown. 
The accord between the curves is within a few percent on-axis 
(triangles) and 6 arcmin off-axis (circles). 
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Fig. 10. Comparison of effective areas of a multilayer-coated 
Wolter-I mirror, as in Fig. |9] but with the source at a finite dis- 
tance (6 = 5 arcmin, D x;102 m). The error bars of the ray-tracing 
outputs are not shown. The accord between the curves is within 
a few percent on-axis (triangles) and 4 arcmin off-axis (circles). 



sist of 200 pairs of Pt/C layers. The layer thickness decreases 
from the coating surface towards the substrate, according to the 
well-known power-law for the d-spacing - i.e., the sum of the 
thicknesses of two adjacent layers - dj = a{b + jy^ (Joensen 
11995 1), with 7=1,2,..., 200 and a, b, c, parameters with values 
depending on the desired reflectivity. In the present example we 
adopted a - 1 15.5 A, b - 0.9, and c = 0.27, for a constant thick- 
ness ratio of Pt to the d-spacing, F = 0.35. The outermost layer 
is Pt. Finally, the surface roughness of the mirror is assumed to 
have the constant value cr = 4 A. 

The results of this test are shown in Fig. |9]for a source at 
infinity, and in Fig. [TO]for a source at a D ^ 102 m distance 
from the mirror The results for a 10** ray tracing for each energy 



value (1 keV steps) are plotted as symbols, whereas the results 
for the application of Eq. ( [26] l are plotted as lines. The statistical 
error of the ray-tracing results is of a few percent at low energies 
and close to 8% at high energies, where the reflectivity is lower 
We note the excellent matching of peak positions and shapes in 
both cases. A mismatch of a few percent can be observed at low 
energies, even at 6 = 6 - 0, which is probably still related to the 
double cone approximation (Eq. (fT2li). 

From these examples, we also note that, even on-axis, the ef- 
fective area is heavily reduced by the source at a 102 m distance 
from the mirror This occurs because 6 = 0.083 deg is close to 
ffo = 0.106 deg, although still smaller, thus the on-axis geomet- 
ric area is only 12% of the one we would have with D - oa (see 
Eq. (|9])). For the case 6 -Q, the effective area also decreases with 
almost everywhere, as predicted by Eq. ( |32] |. The opposite ef- 
fect is observed for a source at a 102 m distance, in agreement 
with Eq. ( [39l ), which predicts an increase in the geometric area 
with 6 when 6 > q'o/2, as in the case that we considered. Finally, 
we also note how the reflectance features are smoothed out when 
the source moves off'-axis, because of the variation in incidence 
angles over the reflecting surfaces. 

6. Final remarks and conclusions 

We have shown how the problem of computing the effective area 
of a Wolter-I mirror with / » L can be reduced to the compu- 
tation of an integral (Eq. ( l25T l, or Eq. ( l26T l as a particular case), 
on the only condition that we are able to compute r,i{a), the mir- 
ror reflectivity as a function of the photon wavelength and the 
incidence angle. This can easily be achieved numerically for a 
source at infinity, as in astronomical cases, or for a source at fi- 
nite distance, as usually done for on-ground calibrations. For the 
ideal case of a constant reflectivity r - 1, we could solve the 
integral and obtain algebraic expressions for the geometric area 
of the mirror. Finally, we presented some examples of the appli- 
cation of the formalism, and checked that its predictions agree 
with those of a detailed ray-tracing routine. 

The analytical approach undoubtedly has several advan- 
tages. In general, the multilayer reflectivity computation is the 
Achilles' heel of a ray-tracing program that is aimed at deter- 
mining the off-axis effective area of a mirror: although concep- 
tually simple, the reflectivity computation consists of summing 
up the contribution of several layers to the reflectivity, so it re- 
quires a significant amount of computation time. A ray-tracing 
routine often requires 10^-^ 10^ photons to return sufficient statis- 
tics for each photon energy, and the reflectivity routine has to be 
called for each of them. Hence, the total computation time can 
reach several hours. If the optic is still to be designed, the entire 
simulation needs to be run several times, adjusting the parame- 
ter values every time, until the optimal solution is reached: the 
optimization process can thereby take several days. 

In contrast, the analytical approach presented in this work is 
completely unaffected by statistical errors. Its accuracy is lim- 
ited only by the applicability of the double cone approximation 
(discussed in Sect. |2]i and by the accuracy of the computation of 
the integral in Eq. ( |25] l. In practice, the reflectivity can be com- 
puted with, say, a 5 arcsec step of the incidence angle, without 
noticeably affecting the reflectivity of the multilayer. This means 
that, even for a very large off-axis angle, e.g., 10 arcmin, we need 
to compute the multilayer reflectivity only ~ 500 times, at most, 
instead of more than 10'* times as required by a ray-tracing. 

For these reasons, the approach presented in this work might 
be extensively used to compute the effective area of grazing- 
incidence Wolter-I astronomical mirrors, on- and off-axis, when- 
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ever applicable. On the other hand, as already stated in Sect. [T] 
it would not be applicable to very short focal lengths, or to sys- 
tems of densely nested mirror shells, obstructing each other in 
the field of view: in this case a ray-tracing is, as of today, the 
only viable computational technique. Clearly, a ray-tracing is al- 
ways necessary to investigate the angular resolution, especially 
whenever mirror deformations are present. 

We note that, if analytical expressions for r^{a) were avail- 
able, it would be possible to solve explicitly the integral in 
Eq. jZSl l in the most general case, and obtain algebraic formu- 
lae for the effective area of a Wolter-1 mirror. In this respect, 
an analytical approach for obtaining r^{a) of a multilayer was 
developed by Kozhevnikov et al. ( 1200 11 1, but a concise expres- 
sion for the reflectivity appears still to be unavailable. However, 
if this is achievable, the computation of the effective area, and 
consequently the optical design and optimization, for any astro- 
nomical X-ray mirrors of sufficiently large /# might be simply 
reduced to the application of a handful of algebraic equations. 

Appendix A: Vignetting for double reflection - 
on-axis source at finite distance 

We derive Eq. (|9l), which returns the geometrical vignetting of 
a double-reflection mirror, caused by the finite distance of the 
source. We consider a radial section of a Wolter-1 mirror (Fig.[T]), 
in the xz plane, with the primary and secondary segment sur- 
faces intersecting at a 2ao angle. We suppose that, as far as the 
effective area is concerned, the mirror profile can be approxi- 
mated by a double cone. If the source were at infinity and on- 
axis (Fig.[T^), the incoming rays would impinge the primary and 
secondary segment at ffo. We now move the X-ray source to a 
finite, although large, distance D. The beam is no longer colli- 
mated, since it has a nonzero half-divergency 6 ~ Ro/D, where 
Rq is the radius at the intersection plane (z = 0) of the mirror 
Because of the smaU cross-section of the mirror as seen by the 
distant source, 6 can be considered as a constant. 

In this configuration, inspection of Fig. [TJ) shows that the 
incidence angle on the primary segment becomes ao + 6, and 
that on the secondary one, ao - 5. We consider a ray striking the 
primary segment at A = {-Rq - Zquq, 0, zo), with zq < L\. The 
angle formed by the ray direction with the optical axis, after the 
first reflection, is 2q'o + 6: the equation of the reflected ray is then 



X + Rq + zpao 
2ao + 6 



■Zq, 



(A.l) 



where, as usual, we approximate the tan functions with the 
respective small arguments. If zo is sufficiently large, the ray 
misses the second reflection (Fig.[T}3). The last reflection on the 
secondary mirror segment occurs if the reflected ray passes by 
the point B - (-Ro + 3a()L2,0, -^2)- Substituting these coordi- 
nates into Eq. (lA.ll ). and solving for zo, we obtain Zq, the maxi- 
mum value of Zo for which we have a double reflection 



ao — 6 

Zo - ^^2, 

ao + 



(A.2) 



regardless of Rq. Then the fraction of the primary segment that 
is effective for double reflection is V = Z^/Li, i.e.. 



V 



L2(aQ - 6) 
Li(aQ+d)' 



(A.3) 



We now consider the entire mirror, obtained by a rotation of the 
profile in Fig. [T] about the optical axis. If the source is still on- 
axis, at a distance D, the vignetting given by Eq. ( |A.3t is easily 



applicable to all sectors of the mirror, and we obtain Eq. (|9]l ex- 
actly. In AppendixiBl we see that this result can easily be gener- 
alized to a source off-axis. 

Appendix B: Incidence angles and vignetting for 
double reflection - detailed calculation for a 
source off-axis at finite distance 

We consider a double cone mirror (Fig. IB. lb . with optical axis 
aligned with z and the intersection plane at z = 0. We define Ro 
to be its radius at z = 0, ao the incidence angle for an on-axis 
source placed at infinity, and Li , L2, lengths of the primary and 
secondary mirror segments. We assume ao to be shallow and 
that the source is at the finite distance D » Li . If the source is 
on-axis, the beam impinges the primary segment with a nearly 
constant half-divergency 6 = Ro/D. We denote with ri, (fx, and 
Zi the radial, azimuthal, and axial coordinates of the impact point 
on the primary segment, and r2,(f2, and Z2 on the secondary. We 
now assume the source to be moved off-axis by an angle 9, and 
choose the direction ip\ - (p2 = in the tilt plane of the source 
(the xz plane, like in Fig.lTJ. Our scope in this appendix is to de- 
termine analytically the incidence angles on the two mirrors, ai 
and a2, and the vignetting for double reflection, V, as a function 
of ao, 6, 6, (fx, more generally than we did in Appendix lAl 

The two conical surfaces are described in polar coordinates 
by the equations 



Zi = with rx > Ro, 



Z2 = 



ao 
r2 - Ro 
3ao 



with r2 < Ro- 



(B.l) 
(B.2) 



The normal vectors to the two segments, directed inwards, are 



•COSQ-qCOS^I ' 

-cosorosin^i 
sin ao 



•±2 



- COS Sao cos ip2 ^ 

- cos 3ao sin ^2 

sin Sffo 



(B.3) 



if the double cone approximation is valid, the normal vectors are 
only a function of the i/j's angles. Ifiht source were on-axis, the 
initial direction of the ray would have the expression 



ha 



' sin5cosi/3i ^ 
sin 5 sin (^1 
-cos (5 



(B.4) 



Since the source is off-axis by Q, we have to tilt by an angle 
Q about the y axis. The application of the rotation matrix returns 
the expression for the initial direction of the oflF-axis photon, 



sin 5 cos if\ cos Q - cos 6 sin Q 

sin (J sin 
- sin 6 cos (^1 sin - cos 6 cos Q 



(B.5) 



If ax is the incidence angle of the first reflection (measured from 
the surface), we can write the scalar product 



cos^^ H- aij = -Hj, (B.6) 

which becomes, after some algebra, 

sinai - cosQ'()Sini5(cos0cos^i^i + sin^i^i) H- 

- cos(q'o -H 5) sin^cos 1^1 -I- cos 5cos ^sin ao- (B.7) 

In the limit of small ao, S, 6, we can approximate the cosines 
with 1 and the sines with their arguments, yielding 



Qfi ^ ao + S — 0COS ifi. 



(B.8) 
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to the focal plane 



Fig.B.l. Double reflection of an off-axis ray (dashed line) on 
a double cone mirror. Also shown are the direction vectors of 
the ray, and the normal vectors to the surface, 
impact points, £j and r^. 



and 112, at the 



that is exactly Eq. ( |20] l. 

To derive the incidence angle on the secondary mirror seg- 
ment, we need to trace the exit direction of each ray, k^, after the 
first reflection. This is obtained from the vector equation 



^1 



nj)nj, 



(B.9) 



because the parallel component to the surface is conserved, 
whilst the normal component reverses its sign in the reflection 
process. By substitution of the Eqs. ( IB.5I ) and ( IB.3I ), one obtains, 
as always in the small angles limit. 



-0 + (6 - 2ai) cos ipi 
(6 - 2ai) sinipi 
-1 



(B.IO) 



A ray reflected at the generic point of the primary segment 
has equation r{t) = r j + tk^, with f > 0. Therefore, the reflected 
ray intersects the secondary segment at a position fulfilling 
the condition r(t) - rj, i.e., ~ Z^i - ^Ki foi" some t. This is 
equivalent to constraining the two vectors to be parallel: 



(B.ll) 



where x denotes a cross product. 

Using both Eqs. ( IB. Il l and ( IB.2b . Eq. ( IB. Ill ) can be devel- 
oped into 3 scalar equations, only 2 of which are mutually inde- 
pendent, i.e.. 



/, 3q'() 

1 H cos ip2 

\ ku 



1 fo 
■2Ro — 3ri 1 H cos^i 



(B.12) 



, 3q'o . 
1 + - — sin (p2 



-i- 2Rq - 3ri 1 -I- - — sin ip\ 



(B.13) 



When 6 - 6 - Q, Eq. ( IB.lOl i infers i = - 2aQ cos (p\ and 
k\y - - 2ao sin tpi . In this case, the solution is barely ip2 = (fii 
and r2 = 4Rq - 3ri. 

We now search for a perturbative solution of these equations, 
of the kind ip2 - (pi + s and r2 = 4-Rq - 3ri + ^, with s ipi 
and ^ «; ri. Substituting these expressions into Eqs. ( IB.121 I 
and ( IB.13I I, neglecting terms in ^e, and using Eq. ( IB. 10b for the 
components of k^ , we obtain a linear system in ^ and e, whose 
solution is 



6(go -gi)(ri - Rq) 
2ao — a\ 
20 sin(/Ji (ri - 7?o) 
(4;?o-3ri)(2Q'o-Qfi)' 



(B.14) 
(B.15) 



For small off'-axis angles, ri - Rq < Lao, 2ao - a\ ~ ao, and 
4Rq - 3ri ^ Rq, s must then be of the order of 20L/Rq or less, as 
expected from the simple argument presented in Sect. 13.21 

We can now derive the incidence angle on the secondary mir- 
ror segment, 02- To this end, we write the scalar product 



cos 



(B.16) 



where is provided by Eq. ( IB.31 l. This would require the com- 
putation of cos((/?i + s) and sm{(fi + s), which would yield a 
complicated expression. Nevertheless, all terms in s are of sec- 
ond order, so they can be neglected in small angles approxima- 
tion. We are then allowed to assume that ^1 ^ (p2 and obtain after 
some passages Q'2 - 2ao - Ci; that is. 



02 — ao — 6 + 0COS ifi 



(B.17) 



so we have found Eq. (ISTT l. This concludes the calculation of the 
incidence angles for an off-axis source. 

Finally, we proceed to compute the vignetting factor for dou- 
ble reflection, in the general case of an off'-axis source. In the fol- 
lowing, since (fi ^ (p2, we suppress the subscript and denote with 
(fi their nearly-common value. From the value of ^ (Eq. ( |B.14| l), 
and using Eq. dB.lVb . we derive ^2, 



r2 



2(ao + a[)Ro - 3airi 



(B.18) 



Reflection on the secondary mirror segment occurs if Z2 > -L2, 
i.e., r2 > Rq - 3aQL2 (see Eq. (IB. 2b ). Substituting Eq. dB.lSb into 
this inequality, using Eq. ( IB. lb . and solving forzi, we derive 



zi((f) < Zo = L2 



aiitfY 



(B.19) 



where we indicate explicitly the dependence of the incidence an- 
gles on the azimuthal angle. So the fraction of the primary mir- 
ror segment efi'ective for double reflection, at the polar angle tp, 
is Zo/Li or 



y(^) = 



L2 aiif) 
Li ai(ip)' 



(B.20) 



The last equation generalizes Eq. (IA.3b to the case of a source 
off-axis: it is Eq. ( l22b . which we used for the computation of the 
off-axis effective area in Sect. 13.21 
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